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Abstract 



We calculate the subthreshold production of antiprotons in the Lorentz- 

>. 
OO , covariant RBUU approach employing a weighted testparticle method to treat 

o 

^— ^ I the antiproton propagation and absorption nonperturbatively. We find that 

O . 

tH- , the antiproton differential cross sections are highly sensitive to the baryon 

a^ 

and antiproton selfenergies in the dense baryonic environment. Adopting 
the baryon scalar and vector selfenergies from the empirical optical potential 
rH , for proton-nucleus elastic scattering and from Dirac-Brueckner calculations 



at higher density p > pQ we examine the differential antiproton spectra as a 



H ' function of the antiproton selfenergy. A detailed comparison with the available 

experimental data for p-nucleus and nucleus-nucleus reactions shows that the 
antiproton feels a moderately attractive mean-field at normal nuclear matter 
density po which is in line with a dispersive potential extracted from the free 
annihilation cross section. 



■^Supported by BMFT, GSI Darmstadt and KFA Jiilich. 



I. INTRODUCTION 

The production of particles at energies below the free nucleon-nucleon threshold ('sub- 
threshold production') constitutes one of the most promising sources of information about 
the properties of nuclear matter at high densities since the particles are produced predomi- 
nantly during the compressed stage at high density |T|-^. Antiproton production at energies 
of a few GeV/u is the most extreme subthreshold production process and has been ob- 
served in proton- nucleus collisions already more than 20 years ago [§|-0]. Experiments at 
the JINR 1^ and at the BEVALAC ^,TU\ have provided, furthermore, first measurements 



of subthreshold antiproton production in nucleus-nucleus collisions. Since then the problem 
was taken up again at KEK [ll|] and GSl |]l2l with new detector setups. Various descriptions 



for these data have been proposed. Based on thermal models it has been suggested that 
the antiproton yield contains large contributions from AA^ ^ p + X, AA — > p + X and 



PP ~^ P^ production mechanisms |13-15|. Other models have attempted to explain these 



data in terms of multiparticle interactions |16 



In a first chance nucleon-nucleon collision model (assuming high momentum tails con- 



sistent with data on backward proton scattering) Shor et al. ^7j succeeded in reproducing 
the antiproton yield for the proton-nucleus case; however, these authors underestimated the 
yield by more than 3 orders of magnitude for nucleus-nucleus collisions. This problem was 



partly resolved by Batko et al. ||18l who performed the first nonequilibrium p-production 
study on the basis of the VUU transport equation. Essential for this success was that in A 
+ A reactions the dominant production channel proceeds via an intermediate nucleon reso- 
nance which allows to store a sizeable amount of energy that can be used in a subsequent 
collision for the production of a pp pair. Later on, these findings were also confirmed by 
Huang et al. [ pJ|l . 

These results have led to the suggestion that the quasi-particle properties of the nucle- 
ons might be important for the p production process which become more significant with 
increasing nuclear density. Schaffner et al. [EO] found in a static thermal relativistic model, 



assuming kinetic and chemical equilibrium, that the p-abundancy might be dramatically 
enhanced when assuming the antiproton selfenergy to be given by charge conjugation of 
the nucleon selfenergy. This leads to strong attractive vector selfenergies for the antipro- 
tons. However, the above concept lacks unitarity between real and imaginary part of the 
p-selfenergy and thus remains questionable. Besides this, even in the a — cj-model the Fock 
terms lead to a suppression of the attractive p-field [^ , such that the production threshold 
is shifted up in energy again as compared to the simple model involving charge conjugation. 
Furthermore, the assumption of thermal and especially chemical equilibrium most likely is 
not fullfilled e.g. in Si + Si collisons around 2 GeV/u ||22|| . 

First preliminary results of a fully relativistic transport calculation for antiproton pro- 
duction including p annihilation as well as the change of the quasi-particle properties in the 
medium have been reported in p^. There it was found that according to the reduced nu- 
cleon mass in the medium the threshold for p-production is shifted to lower energy and the 
antiproton cross section prior to annihilation becomes enhanced for Si + Si at 2.1 GeV/u 
by approximately a factor 70 as compared to a relativistic cascade calculation where no 
in-medium effects are incorporated. However, all these transport calculations |T^,|TP|,^ suf- 



fered from an approximate geometrical treatment of the very strong annihilation channel 
and a neglect of the momentum dependence of the baryon selfenergies. 

In our present work we therefore analyze the production of antiprotons in the framework 
of the relativistic transport theory (RBUU) where antiprotons are propagated explicitly in 
the respective time dependent potentials and their annihilation is calculated nonperturba- 
tively by means of individual rate equations. 

First results have been published in [^ . Results of a very similar investigation by Li et 
al. |2^ have recently become available and we compare frequently with the outcome of this 
study. 

In this respect we first discuss two simple models for the p selfenergy in section 2 and 
present the general transport equations for nucleons and antinucleons in section 3. The 
numerical implementation of p production is presented in section 4 as well as a detailed 



analysis of the production process with respect to the selfenergies employed, the systematics 
with respect to projectile and target masses and the individual baryonic production channels. 
In section 5 we present the treatment of p propagation and annihilation and discuss the 
geometrical aspects of p absorption. The explicit comparison of our calculations with the 
available experimental data ioi p + A and A + A reactions is performed in section 6, while 
a summary on the p selfenergies concludes the paper in section 7. 

II. MODELS FOR THE ANTIPROTON SELFENERGY 

The dynamics and the properties of particles in a many-body system strongly depend on 
their mutual interactions with the surrounding particles and are reflected in their selfenergies. 
While the real part of the selfenergy describes the change of the particle momenta in the 
medium, all inelastic reaction channels as well as elastic scattering processes are accounted 
for by the imaginary part of the selfenergy S(p^,ps;J^)- ^^ ^^^ selfenergy S is an analytic 
function of the particle energy e = Pq its imaginary and real part are related up to a constant 
by the dispersion relations: 

/,.iEwi^-ipr^?ii^*'. (2) 

TT JO e — e 

Since we assume particles and antiparticles to be different particle species the integration in 
the Principal-value integrals ([l|) and (^) extends only over positive energies. 

In the relevant energy regime (1-5 GeV/u) of heavy-ion collisions and proton- nucleus 
collisions the selfenergies of baryons and antiprotons are known to be quite different. The 
p-selfenergy is dominated by the antiproton annihilation with baryons {p + B ^ X) as 
shown by low energy antiproton-nucleon scattering experiments at LEAR p6| , p7| . Since 
the antiproton production probability at subthreshold energies is extremly low the number 
of antiprotons produced in heavy-ion collisions and proton-nucleus collisions is negligible 
compared to the number of baryons involved in the reaction. Thus contributions from 



elastic p-p-scattering to the antiproton selfenergy can be neglected as well as contributions 
of p-annihilation channels to the baryon selfenergies . 

A. A dispersive model 

In the following we present a simple model for the real part of the p-selfenergy using 
the dispersion relation (|I]). In the low density limit the imaginary part of the selfenergy 
(neglecting elastic scattering) is given by the integral : 

2 im s,(x, n,-) = -^ j d'YiB ^^(nj^, n^)/(^> n^) (3) 

where iy(n|^, II^,) is the probability of an antiproton with momentum fl^^ to annihilate 
in a collision with a baryon of effective mass m* and momentum 11^. /(x, IIb) is the 
baryon phase-space distribution function at space-time x whereas the factor 4 in formula 
(0) arises from the summation over the spin- and isospin-degrees of freedom. We simplify 
this expression by considering an antiproton moving in infinite nuclear matter of density po 
(= Q.ni/ fm?) at rest, i.e. n^ = 0. Both the antiprotons and the baryons now are supposed 
to have effective masses equal to the restmass of the nucleon. This leads to the following 
analytic expression for the ^-annihilation probability 



wm, n^) = ^ ', ■ a,,.(m) (4) 



where aabs is the total p-annihilation cross section. Due to the above assumptions the integral 
over the baryon phase-space distribution function and the summation over the spin- and 
isospin-degrees of freedom in eq. (0) can be replaced by a multiphcation of the annihilation 
probability with po- The expression for the imaginary part of the antiproton selfenergy as a 
function of the antiproton energy then reads 



1 mo 2 — 777,2 

Im s,(n°) = --)Ljl_^ . ^^,^(nO) p, (5) 

and the real part of S at density po can be calculated from eq. (|^) by employing the 
dispersion relation (|^). We expect this expression to give a reasonable description of ImT^p 
at least for densities up to po where possible medium corrections to Cats are still small. 



In order to estimate ReT,{e) at po we adopt the parametrization from |T3| for aabs, i-e. 



(Tabs = (To— 7 ^2 ° A2 + ^ (^) 

with the constants (Tq = 120 mb, A = 50 MeV, B = 0.6 and sq = 4m^. Fig. 1 displays 
this parametrization in comparison to the experimental data for the free cross section taken 
from |^8| as a function of s — Sq where s is the invariant energy squared. The result of this 



calculation for i?eS(e) is displayed in fig. 2 as a function of the p-kinetic energy. In order 
to investigate the sensitivity of this result to the parametrization of the absorption cross 
section we used in addition different parameters B (cf. fig. 2) to vary the overall size of aabs- 
For an antiproton at rest we find values for the real part of the selfenergy of —175 MeV up 
to —100 MeV. While we observe a steep rise for small kinetic energies this flattens out for 
higher energies and reaches an asymptotic value of — 10 MeV. Furthermore, the variation 
in the cross section mainly affects the low energy regime. The p-selfenergies obtained in 
the framework of this simple model are well in line with the potential analysis for p+ A 
reactions by Janouin et al |^ . 



B. Mean-field models 

We now turn to the treatment of the antiproton selfenergies in mean-field models where 
the imaginary part of the selfenergies is not included at all thus inherently violating the 
dispersion relations (|I|, 0). 



As an example we consider the mean- field approximation of the familiar a — cj-model pO| 
where only momentum independent vector U^{x) and scalar parts Us{x) of the selfenergies 
are taken into account. The equation of motion for a fermion spinor then reads 

{r [td^ - U;{x)) -{m + W {x))] ^ (x) = (7) 

while the equation of motion for antiparticles is obtained by applying the charge conjugation 
operator to this equation. As a consequence the scalar part of the selfenergy is the same 
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for particles and antiparticles while the vector part of the antiparticle potential changes the 
sign 

^c.(^) =-^;(^) (8) 

U^ (x) = W (x) . (9) 

The combination of both parts of the antiparticle selfenergy to the Schr~'dinger-equivalent 
potential in the non-relativistic limit leads to a strongly attractive potential for the antipro- 
tons when using the original coupling constants from pO]. The value of the Schrodinger- 



equivalent potential at p = po is approximately —700 MeV for e = and becomes even more 
attractive with increasing kinetic energy of the antiproton (—1300 MeV for e = 2000 MeV). 
A comparison of the results for the Schr^'dinger-equivalent potential obtained in the 
framework of the a — tu model and the real part of the p-selfenergy resulting from the 
dispersive approach shows strong differences not only in the absolute values but also in their 
opposite behaviour as a function of the kinetic energy. The fact that the real part of the 
antiproton selfenergy cannot be described consistently within different models urges us to 
treat it as a free parameter in our transport model. We regard the determination of this 
parameter in comparison with the experimental data as the central goal of our work. 

III. THE RBUU TRANSPORT APPROACH 

In this section we give a brief description of the RBUU-model. First we summarize the 
relevant equations determining the dynamics of baryons and then describe the implementa- 
tion of antiprotons in our approach. 

A. Baryon dynamics 

Since the covariant BUU approach has been extensively discussed in the reviews pl|j3^ 



and in |]53| we only recall the basic equations and the corresponding quasi-particle properties 



that are relevant for a proper understanding of the results reported in this study. The 



relativistic BUU (RBUU) equation with momentum-dependent mean-fields or selfenergies 
is given by (for details see refs. I31|-p3|) 



{[n, - U^icF^Un + m^idf^UsM + i-^Ad^.Un + m\d^^Us)]d^}fix,p) = Ui, (10) 

wfiere f{x,p) is tlie Lorentz covariant pliase-space distribution function, Icoii is a collision 
term (cf. eq. (^)), and Us and t/^ are the scalar- and the vector- selfenergies. The effective 
mass M* and the kinetic momentum 11^ are defined in terms of the fields by 

^fM{x,p) =p^,-U^{x,p) (11) 

^^^{XjP) = m + Us{x,p) (12) 

(13) 

while the quasi-particle mass-shell constraint is obtained as 

V{x,p)f{x,p) = (14) 

with the pseudo potential 

V{x,p) = ^{Il\x,p) - m*\x,p)}. (15) 

The above equation implies that the phase-space distribution function f{x,p) is nonvanish- 
ing only on the hypersurface in phase-space defined by V{x,p) = 0. 

In order to implement proper selfenergies for the nucleons in line with elastic proton- 
nucleus scattering we follow ref. [Q and separate the mean-fields into a local part and an 
explicit momentum-dependent part, i.e. 

[/,(x,p) = f/f(x) + f/f^(x,j9) (16) 

t/^(x,p) = [/f(x) + f/;^^(x,p), (17) 

where the local mean-fields are determined by the usual Hartree equation 



U"{x) = -gsfrnix) 

U;^ (x) = g^uj! (x) (18) 



with 



mlanix) + 5^o-|(x) + Cscr|^(a;) = QsPsix) 

^l^ja^) = QvUi^)- (19) 

In the above equations the scalar density Ps{x) and the baryon current j^{x) are given in 
terms of momentum integrals over the phase-space distribution function (cf. [p^ -p5[]). The 
potentials (|r^) and (|TB|) correct the unphysical strong repulsion in the p + A and A + A 
potentials at high erngergies, obtained in the original Walecka model. Refs. p3i|34| give 
parametrizations which describe both the experimental data and the density-dependence 
obtained in Brueckner-calculations very well. 

Now we turn to the discussion of the collision term Icou describing the baryon-baryon 
collisions (cf. |[1|. |31|J36|] ): 

Icon = -, — 7T / rf^Hi rf^n rf^n'^ — — -i- 

w{w, n^ I n''^, n'/')5^(p^ + p^ - p'^ - p'f ) 

(/(x, n')/(x, n;)(i - /(x, n))(i - /(x, no) 

-/(x, n)/(x, ni)(i - /(x, n'))(i - /(x, n;))). (20) 

This collision integral describes the change in the phase-space distribution function /(x, 11) 
due to the collision of two baryons with effective masses m*, m\ and momenta 11'^, J\^^ 
respectively. The two baryons in the final state of the reaction with masses m'* and m'^^* are 
labeled by their momenta 11' and 11']^. The 5— function guarantees energy and momentum 
conservation in the individual collision while 14^(11'^, W^ \ 11'^, II'/^) denotes the transition 
probability for this reaction which can be expressed in the CMS of the colliding particles 
by the product of the relative velocity of both colliding particles and the differential cross 
section for the reaction, 
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W{W, n^ I ^'^ n^) =1 ^ + ^ I ^ In.+n^'^n'.+n- • (21) 



For the NN cross section we use the parametrization given by Cugnon [^ and therefore 
neglect the possible in-medium corrections. However, the corrections due to Pauli-blocking 
are built in via the factors 1 — /. We explicitly account for the following baryonic channels: 

and propagate the A's with the same selfenergies as the nucleons. For the parametrizations 
of the cross sections of reactions including A's see [|l[. 



B. Antiprotons in RBUU 

The phase-space distribution function for the antiprotons fp{x,Pp) is assumed to follow 
an equation of motion equivalent to eq. (0), however, with scalar and vector potentials of 
different strength, i.e. 

U!{x) = -g^anix) 

Kix) = glu^'^ix),. (22) 

Thus the effective mass and the effective momentum of antiprotons are given by 

ml = m + UP{x) =m- glan^x) (23) 

K=r^,- Ui{x)= Pi - glu^ix),. (24) 

According to the arguments given in section |T| the coupling constants (^f and g^ are treated as 
free parameters and will be determined in a comparison of our calculations to experimental 
data (cf. section [V| ). Since at subthreshold energies antiprotons can be produced only 
in a very limited kinematical range we can work here with scalar and vector parts of the 
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antiproton selfenergies that are not explicitly momentuin-dependent. We then obtain the 
following equation of motion for the antiproton phase-space distribution function 



' 5,^ + K^^I. + ^;(5>;-))5{;, 



n%- 



U{x,U,) = Il,,{x,Up) (25) 



with 



Ff = d^UP'' - -d'^Ui^ (26) 

and the mass-shell constraint 

(n|-mf)/(x,n,-) = o. (27) 

The collision term I^^n (r.h.s. of eq. (^)) includes i) a term Ili^sti describing the elas- 
tic antiproton scattering, ii) a term Ip^^di describing the antiproton production, and iii) a 
term /^^^^ responsible for in-medium p-absorption. /f^^g^ describes elastic baryon-antiproton 
scattering as well as elastic antiproton-antiproton scattering. While this part of the colli- 
sion integral I^^^n can be formulated similarly to the collision term describing bayon-baryon 
scattering (^) the other terms represent extensions to this integral. 

1. Collision term for p-production 
The basis for the description of the antiproton production is the reaction 

B + B^p + p + N + N=l + 2^p + ?> + A + b, (28) 

for which the corresponding covariant collision integral reads 



* ^^* »^* ^^*- »^* .v^* 



(27r)iiy 12 3 4 5 -^ononon^ngno- 
PFin^ n^ I n^, ui n^, w^)6\p'i + ^^ _ ^^ _ ^^ _ ^m _ ^,^ 

{fix, ni)/(x, n2)(i - fix, n3))(i - fix, n4))(i - fix, u,))} . (29) 
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We have omitted the Pauh-blocking factor for the antiproton in the final state because the 
number of antiprotons created during a heavy-ion coUision in the subthreshold energy regime 
is negligible. For the same reason we neglect the effects of reaction ( pH]) on the phase-space 
distribution function of the baryons. 

2. Collision term for p-absorption 

In the collision integral /^^^ we do not treat all possible annihilation reactions separately 
but sum up all channels in the inclusive annihilation reaction 

B + p-*X, (30) 

where X denotes all possible final states (essentially pions) of the baryon-antiproton annihi- 
lation. The corresponding energy and momentum conservation reads 

p'b+p'p=p'x. (31) 

where p^ and pp denote the 4-momenta of the baryon and the antiproton, respectively, and 
Px stands for the sum of the 4-momenta of all particles in the final state of the annihilation 



reaction. Reaction (pOD then leads to the collision term 



P rr. TT_^ — / ^3tt. ^4tt 1 P 



Il,X^.Iip) = --— ^ / d^Ii^d'^Yix-^ 



ly (n^ n^ I Wx¥\p^i +p'p- p'x)fi^, ni)/(a;, n^-) (32) 



with 1^(11]', lip I n^) denoting the transition probability for the reaction (|30|) . Integrating 
(^) over d^IIx implies in addition to the integration over all final momentum states of a 
particular reaction a summation over all possible annihilation channels 

/L(a;,n,) = -^_^|d3n^^i^IV(n^n^)/(a;,ni)/(a;,n,), (33) 

where 1^(11]*, 11^) denotes the probability of an antiproton with effective momentum 11^ 
annihilating with a baryon with effective momentum Hi. 
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IV. ANTIPROTON PRODUCTION IN RBUU 

In this section we describe the numerical treatment of antiproton production in our model 
and present a systematic analysis of the p-production mechanism in heavy-ion collisions . 

A. Numerical implementation 

Since the production probability for antiprotons is very small the average time evolution 
of the nucleus-nucleus collision is not affected and it is justified to treat the p-production 
perturbatively P,pT|. In this approach the antiproton invariant differential multiplicity is 
obtained by summing incoherently over all baryon-baryon collisions and integrating over all 
residual degrees of freedom. Assuming the antiproton production to take place via reactions 
of the type 

B + B^p + p + N + N = l + 2^p + 3 + 4: + 5 (34) 

(B stands for either nucleon or A) the invariant multiplicity as a function of the impact 
parameter can be written as 



d'^p B'k^oiiJ Vbb(v^) ^(PYi'^d?n',dm',d^w^ 

(1 - /(x, n'3))(i - /(x, n;))(i - /(x, n's)), (35) 

where the quantities IIj (i = 1,..,5) denote the in-medium momenta of the participating 
baryons. lip and Ep stand for the p effective momentum and energy while s = {Hi + 
112)^ is the squared invariant energy available in the corresponding baryon-baryon collision. 
An integration over the impact parameter then yields the Lorentz invariant differential 
production cross section 

We assume the elementary antiproton production cross section crBB^p+xiy/s) in all baryon- 
baryon channels to be equal to <Jpp-^p+x{^/s) and employ a parametrization of the free cross 
section given by |jl8| (cf. fig. 3; solid line): 
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<yBB^p+x = 0.01(v^ - v^)^-«^^[m6], (37) 

with yso = 4m7v and mAr = 0.9383 GeV/c^. The dashed and the dotted hne in fig. 3 
represent extreme alternative parametrizations that will also be used below. 

While in free space the threshold for the elementary production reaction is obviously 4 
times the nucleon restmass in the medium one additionally has to take into account the 
selfenergies of all participating particles. The conservation of energy and momentum has to 
be guaranteed, i.e. 

Pt+P'i=p't + Pt + P't+Pl, (38) 

which in terms of effective momenta and selfenergies (cf. eqs. ( pTj) and (12^) ) yields 



n^ + u^{\ pi I, x) + n^ + u^{\ P2 1, x) = n^ + f/,^(| ^3 1, x) + n^ + 

ui^il P4 I, x) + n^ + ui^ii p5 Ix) + ^ + U^^ix). (39) 

With the abbreviation 

A^ = Ui^il p3 |,x) + U>;i\ p4 |,x) + Ui^il p, \,x) + Urix) 

-U,^{\pi\,x)-U{:{\p2\,x) (40) 

we obtain in shorthand form 

n^ + n^ = n^ + n^ + n^ + n^ + a^. (41) 

When evaluating the threshold for the reaction in terms of effective momenta, i.e. 

tip + lis + n4 + lis = 0, (42) 

we in general encounter a nonvanishing sum of the vector selfenergies of all participating 
particles. This leads to the following expression for the threshold 



^ = J (mt + m^ + ml + mt, + A^)'^ - IS?. (43) 
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In the a — cu-inodel, where the p-selfenergies result from the corresponding baryon selfen- 
ergies by applying charge conjugation, the quantity A'^ vanishes and eq. ( ^31) reduces to 
an expression equivalent to the free threshold with the restmasses replaced by the effective 
masses of the participating particles. 

In order to derive an expression for the differential p-multiplicity we assume, as in refs. 
ITgHTsI, that the differential elementary antiproton production cross section is proportional 



to the phase-space available for the final state: 

hah^h^hp - aBB~.NNN+pW-V 



5^(n^ + n^-n|^-n^-n^-n^-A^). (44) 

Here, the 5-function guarantees the energy and momentum conservation and ^/s is the 
invariant energy available for the quasi-particles in the initial state. R4,{^/s) is the 4-body 



phase-space integral |38]; it has been included to ensure that the differential cross section is 



normalized to the total cross section. 

B. Sensitivity to the elementary cross section 

Now we turn first to the analysis of the antiproton production mechanism itself and 
consequently neglect all in-medium propagation and absorption effects. For the sake of nu- 
merical simplicity we calculate here the p-production using the baryon-selfenergies obtained 



from the nonlinear a — cu-model |39| where the p- self energies are determined by charge con- 
jugation from those of the baryons. 

Since we are dealing with subthreshold particle production the antiproton cross section 
will depend strongly on the behaviour of the elementary cross section ([37| ) close to threshold. 

Since there are no experimental data available for ^/s ^/so < 1 GeV (cf. fig. 3) we have 

to rely on an extrapolation of the data to threshold. In order to investigate the dependence 
of our results on the parametrization of the cross section we display in fig. 4 the invariant 
differential production probability for the reaction Si + Si at 2.1 GeV/u for different extreme 
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parametrizations corresponding to the dashed and dotted hnes in fig. 3. This uncertainty 
of up to one order of magnitude has to be kept in mind when drawing any conclusions from 
the comparision of our work with the experimental data. For our further analysis we will 
adopt the parametrization (pTf). 

Due to these uncertainties it is justified to neglect the Pauli-blo eking effects in the cal- 
culation of antiproton production which turn out to reduce the production probability by 
about 10 to 20% as shown in fig. 5 for the special case of a central Au + Au collision at 2.1 
GeV/u. The effects of the Pauli-blocking increase with decreasing beam energy and become 
less important for lighter systems as e.g. A^e + Ne. 

C. Probing the high density phase ^vith p 

In order to demonstrate the sensitivity of the antiproton production to the equation of 
state (EOS) for nuclear matter we display in fig. 6 the invariant differential p-production 
probability for the reactions Au + Au and Si + Si at 2.1 GeV/u calculated for different 
parametrizations of the EOS. The parametersets NLl and NL3 (cf. ref. ||3^) employ the 
same incompressibility K = 400 MeV, but differ in the effective nucleon mass at p = po 
{m* /m = 0.83 for NLl; m*/m = 0.7 for NL3), while NLl and NL2 employ the same effective 
mass, but differ in their incompressibility {K = 200 MeV for NL2). For both systems we 
observe a production probability which is enhanced by roughly one order of magnitude for 
NL3 as compared to that from NLl and NL2. The reason for this behaviour is obviously 
the low lying threshold which goes along with the small effective nucleon mass when using 
the prametrization NL3. 

While there is hardly any sensitivity to the incompressibility (NLl versus NL2) in the 
light system, one can observe an increase of the production probability by a factor of 3 with 
decreasing incompressibility in the heavy system. The reason for this sensitivity is given in 
fig. 7 where we show the differential p-multiplicity dN/dp as a function of the baryon density 
p/po at the individual antiproton production point for Si + Si (dashed line; multiplied by 
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a factor of 50) and Au + Au at 2.1 GeV/u. It is clearly seen that almost all antiprotons 
are produced at p > 3po in the heavy system while the majority of antiprotons is created 
between 2po and 3po in the light system. This difference is easily understood in terms of the 
pile-up of density in the heavy system. In the case of Au + Au reactions high densities are 
obtained for parametersets that describe a rather soft EOS HO]. Thus compared to NLl the 
parameterset NL2 leads to a larger density in the reaction zone resulting also in a higher 
antiproton production cross section. 



D. Variation with the impact parameter 

Since central collisions lead only to minor contibutions to the inclusive cross section in 
heavy-ion collisions it is of interest to study the p-multiplicity as a function of the impact 
parameter b. For this purpose we show in fig. 8 the calculated differential antiproton 
multiplicity for Si + Si at 2.1 GeV/u multiplied by 27Tb as a function of b (solid line). This 
quantity can be approximately fitted by (dotted hne) 

da{b) 



db 



2TTbPib = 0) S{b), 



(45) 



where S(b) represents the geometrical overlap of two nuclei with mass number Ai and A2 
as a function of the impact parameter. Assuming Ai > A2 the overlap function S(b) reads 
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{¥ + Ri- Ri) + {¥ + Ri- Ri) J2{Rl + Rl) b^-b^- [Rj - R^)^ 



for R2<b<Ri 



and 



S{b) = 1 
for b < R2 < Ri 



(46) 



with Ri = aAl and R2 = aA2 . In case of symmetric systems this reduces to 
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A more systematic analysis with respect to Ai and A2 shows that an optimal fit to the 
impact parameter dependence is obtained for the radius parameter a = 1.1 fm. 



E. Baryonic decomposition 



It is well known from transport calculations |]18| , ^ that heavy mesons and antiprotons 
are dominantly produced via multiple baryon-baryon collisions where especially nucleon- 
resonance channels play an important role. To gain further insight in the mechanism of 
p-production we analyse the relative contribution from various channels in more detail by 
counting the number of collisions A^^i + A^2 that two baryons have undergone before producing 
an antiproton in a mutual collision and determine the corresponding production probability 
as a function of A^i + N2 for the channels A^A^, A^A and AA separately. The resulting 
distributions are displayed in fig. 9 for the reactions Si + Si, Ca + Ca and Au + Au at 2.0 
and 2.5 GeV/u. The full histograms denote the channel A^A^ —^p + X, the open histograms 
denote the channel A^A -^ p+X and the hatched histograms the A A -^ P+X, respectively. 
Table | shows the relative contributions integrated over A^i + A^2- From table 1 we can extract 
a general trend: The A^A^-channel is obviously of minor importance. In heavier systems the 
contributions from resonances to the p-production are generally higher as compared to lighter 
systems. This correlates with the higher density pile-up in heavy systems. The resonances 
become, furthermore, more important with decreasing beam energy. 

In order to study this phenomenon in more detail we display in fig. 10 the relative 
contributions to the p-yield for Ca + Ca as a function of the beam energy. The contribution 
from the A^A^-channel is below the 10% level for Eheam ^ 2.56*6^/^. While for higher 
energies the AA^-channel dominates the production of antiprotons the AA- and the AA^- 
channel are approximately equally important in the energy regime between 2.4 GeV/u and 
2.0 GeV/u. Below 2.0 GeV/u the AA-channel is clearly the dominant one. 



From the analysis of this section we can conclude that, as far as antiproton production 
is concerned, the nucleon-resonances act as energy reservoirs that can release their energy 
in baryon-baryon collisions to allow for p-production. That the resonances are populated 
in sufficiently large numbers is a consequence of the formation of resonance matter ^3[ in 
relativistic heavy-ion collisions. 

F. Sensitivity to in-medium p-properties 

The threshold for the elementary p-production reaction according to (|43|) depends on 
the in-medium properties of the antiprotons . In order to investigate this dependence we 
calculate the p-production probability for Ni + Ni at 1.85 GeV/u (central collision) for 
a variety of different antiproton selfenergies : In the calculation shown in the upper part 
of fig. 11 the p-production probability was determined using a fixed vector part of the 
antiproton selfenergy {Uq = —98 MeV at p = po) while varying the scalar part f/f (see 
figure caption). The results displayed in the bottom part of fig. 11 were obtained by varying 
the vector part while keeping the scalar part f/f fixed at —159 MeV at p = po. The same 
analysis performed for other systems than Ni + Ni lead to quite similar results. Thus we 
conclude that the antiproton production cross section exhibits an extreme sensitivity to the 
in-medium p-properties. This sensitivity will be used below to obtain approximate values 
for the p selfenergy in comparison to experimental data. 

V. IN-MEDIUM P PROPAGATION AND ABSORPTION 

The antiprotons produced in individual baryon-baryon collisions can be annihilated while 
propagating out of the dense nuclear medium. Since we expect the p-annhilation to dominate 
the antiproton spectra a proper treatment of this reaction channel is crucial for a comparison 
with experimental data. In order to be able to treat the p-absorption and propagation 
sufficiently well we employ a new numerical method which allows to calculate the propagation 
as well as the antiproton absorption nonperturbatively. We neglect the elastic scattering 
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of antiprotons with the surrounding baryons; according to ref. p5| this elastic scattering 
only causes a flattening of the antiproton momentum spectra. The numerical details of 
our method are presented in Appendix A. In this section we concentrate on the effects of 
in-medium absorption and propagation of the antiprotons . 

During a heavy-ion collision or a proton-nucleus reaction the antiprotons are created 
with effective masses m| and momenta Up. These particle properties change according to 
the surrounding medium as the antiprotons leave the reaction zone. At the end one observes 
free antiprotons with kanonical momenta pp. Since the properties of the antiprotons are 
determined by their optical potential one expects signiflcant effects of the p-propagation on 
the p-spectra. We analyse these effects in flgs. 12 and 13. Fig. 12 shows the variations in 
the p-spectra due to different antiproton selfenergies for a central Si + Si collision at 2.1 
GeV/u at Qiab = 0°. The production probability is identical for all cases. The calculations 
shown in the top part of this flgure were performed for f/f = and different vector parts 
(see flg. 12). The bottom part of this flgure shows the result of calculations for f/^ = and 
varying scalar parts of the antiproton selfenergy . Comparing both flgures we conclude, as 
for the p-production itself, that the effects of the scalar and vector part of the selfenergy 
on the p-spectra are qualitatively and quantitatively similar. Attractive potentials result 
in a reduction of the p-momenta as the antiprotons leave the reaction zone while repulsive 
potentials lead to enhanced p-momenta. 

Fig. 13 displays the same analysis for p + Cu at 4.0 GeV at Qiah = 0°. Again, the 
production was calculated identically for all three cases. Comparing the p-spectrum obtained 
for free propagation {Ugep = 0) with the one obtained using a weakly attractive optical 
potential {Usep = —90 MeV) one observes a slight shift of the spectrum towards lower 
momenta. This shift is easily explained by the loss of kinetic energy that occurs as the 
effective mass nip increases to the value of the restmass when the antiprotons leave the 
reaction zone. Using a strongly attractive optical potential {Usep = —560 MeV) one reahzes 
in addition to this shift a reduction of the spectrum due to the propagation effects. Contrary 
to heavy-ion collisions the target-nuclei in proton-nucleus collisions are not destroyed during 
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those tiinescales when the antiprotons move into the continuum. This means that the target 
nuclei build a potential barrier which the antiprotons have to overcome when leaving the 
nuclei. The low energy antiprotons cannot overcome this potential wall and get trapped in 
the nucleus which leads to the observed reduction in the p-spectrum. 

Figs. 12 and 13 clearly show that the antiproton in-medium propagation has an impact 
on the differential p-spectra. Therefore it is important to apply the complete mean-field 
dynamics of the RBUU model when describing the propagation of antiprotons . 

Now we turn to the in-medium antiproton absorption {p + B -^ X). Due to the fact that 
there is no information available on the in-medium cross section we have to rely on the free 
cross section in the parametrization from |T^(6). In our transport calculation we replace 
the free s — Sq with that determined from the in-medium properties of the colliding particles, 
i.e. s = {Up + n^)^ — (lip + n^)^ is the squared invariant energy available in the elementary 
antiproton-baryon collision while Sq = {rrip + m*^)'^ denotes the squared sum of the effective 
masses of the colliding particles. The dashed line in fig. 1 represents a cutoff at 100 mb which 
is introduced to simulate possible in-medium screening effects. The value of 100 mb for the 
cutoff is in line with the annihilation radius for proton-antiproton annihilation derived in a 
optical model calculation [^. Fig. 14 shows the number of antiproton-baryon collisions as 
a function of the corresponding value of aabs for a central reaction Si + Si at 2.1 GeV/u. 
Since only a small fraction of all events lies in the region of 100 — 192 mb our results do not 
depend significantly on the cutoff. Similar studies for proton-nucleus collisions show even 
less sensitivity to this cutoff since the antiprotons move with higher momenta with respect 
to the target nucleus. 

In order to study the p-absorption effects on the antiproton spectra from heavy-ion 
collisions we have performed calculations using different fixed absorption cross sections in 
comparison to the dynamically determined cross section (^. Fig. 15 shows a full calculation 
for the system Si + Si at 2.1 GeV/u at 6cm = 0°. The top line denotes the calculated 
p-probability without absorption using the baryon selfenergies obtained from the nonlinear 
a — cu— model with p- selfenergies obtained via charge conjugation from the baryon selfener- 
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gies. The lower lines represent calculations with different constant absorption cross sections 
while the solid line results from a calculation using the parametrization (|^). This figure 
clearly reffects the dominant role of the in-medium absorption of antiprotons which leads 
to a reduction in j)-probability by approximately two orders of magnitude; these absorption 
rates agree with those determined by Li et al. |^. Furthermore, the spectral distribution 



obtained when using the parametrization (^ is - except for small deviations at low antipro- 
ton momenta - almost identical to the distribution obtained for a constant cross section of 70 
mb. This value corresponds approximately to the mean value of the distribution displayed 
in fig. 14. 

In this context we note that in proton-nucleus reactions the corresponding mean value 
of the elementary p-absortion cross section is aproximately 50 mb which leads to absorption 
factors of 10 — 30 depending on the size of the target nucleus and the beam energy. The 
difference in antiproton absorption between heavy-ion collisions and proton-nucleus reactions 
is mainly due to the different kinematics in both cases: During a heavy-ion collision most 
of the antiprotons are created in the CMS of the collision which implies that the relative 
velocity of the antiprotons and the surrounding baryonic medium is low. This leads to a 
small s — So for the annihilation reaction and results in high values of the absorption cross 
section (see fig. 1). In contrast to this situation the antiprotons produced in proton- nucleus 
collisions move with momenta of around 1 GeV/c through the baryonic medium which 
leads to smaller annihilation rates. Another difference in both types of reactions is that 
the p-absorption in proton-nucleus collisions takes place at a maximum density po while the 
antiprotons produced in heavy-ion collisions experience higher densities (cf. fig. 7) which 
again leads to an increase in the absorption probability. 

Finally we study the impact parameter dependence of the antiproton absorption. For this 
purpose we have to investigate a physical quantity that does not show the impact parameter 
dependence of the antiproton production mechanism. The ratio R of the differential p- 
probability calculated including absorption and the differential production probability meets 
this requirement. In Fig. 16 we display R for Ni + Ni at 1.85 GeV/u, Ne + Ni at 2.0 GeV/u 
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and Si + Si at 2.1 GeV/u as a function of b/bmax where b denotes the impact parameter and 
bmax is the sum of the corresponding radii of target and projectile. The ratio R is evaluated 
for zero antiproton momentum in the CMS of the corresponding heavy-ion collision . The 
dots represent the values of R calculated in the RBUU model and the solid lines correspond 
to fits to these numerical data using the function 

R{b/b^ax) = A-{l + e''''''^^)-\ (48) 

The parameters A, B and C are given in table ||. Except for peripheral collisions {b/bmax > 
0.6) the ratio R is determined by the constant A for all systems. This means that the 
p-production spectra for b/bmax < 0.6 are reduced due to the absorption mechanism by a 
constant factor determined by the size of the system. The parameters B and C are similar 
for all systems considered. This leads to the conclusion that the absorption mechanism is 
similar for all systems and can be understood by means of simple geometrical considerations. 

VI. COMPARISON WITH EXPERIMENTAL DATA 



After the systematical analysis of the in-medium antiproton production and absorption 
mechanism we now turn to the comparison of our calculations with the most recent data 
from KEK and GSI. In order to perform this comparison we describe the p-selfenergies on 
the basis of the a — a;— model with free coupling parameters g'^ and g'^. The comparison 
with the experimental data will allow to approximately determine these parameters and then 
provide first information on the antiproton potential in the dense medium. The quasi-particle 
properties, i.e. the nucleon selfenergies Us{x,p),Uf^{x,p), of the baryons participating in 
the p-production reaction are taken from refs. p3|-[35[l. Us{x, p) and U^{x, p) are fixed to 
reproduce the saturation properties of nuclear matter, the empirical proton-nucleus optical 
potential as well as the density dependence of Us and f/^ from Dirac-Brueckner theory [^ . 



For orientation the actual values for Us{p) and the zero'th component of the vector field 
Uq{jp) are displayed in fig. 17 for po(~ 0.17/m^^), 2po, and 3po- 
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Before we present the results of the calculations using these expressions for the baryon 
selfenergies we give a brief description of the numerical method used to implement these 
selfenergies into the antiproton production process. Since we deal with explicit momen- 
tum dependent selfenergies {Us{x, p), U^{x, p)) for the baryons their effective momenta and 
masses are given by 

^=P^-Ui^i\p,lx) (49) 

m* = m + Us{\pj \,x) (50) 

J = 1, ■■■, 5, 

while the corresponding quantities for the antiprotons read 

^ = p'p- urix) (51) 

m; = m + f/f(x). (52) 

The elementary p-production events occuring at location x are evaluated in the correspond- 
ing local rest frame (LRF) of the nuclear matter (j^ = (p, 0, 0, 0)). In this frame the spatial 
components of the vector selfenergies vanish by definition. This implies that the vector 
components of the effective and the canonical momenta of all particles are the same (]49| , 



Slf ). The spatial components of the quantity A'^ defined in eq. (|40| ) also vanish. Energy and 



momentum conservation then yields 

n° -f n° = n° + n° + n° + nj- + a°, (53) 

Pl+P2=P3+Pi+P5+ Pp- (54) 

In order to calculate the p-production probability we first have to evaluate the threshold 
y/s^ in the LRF. In the CMS of the particles in the final state the threshold is obtained 
for the kinematical situation with all particles at rest. Thus in the LRF these particles 
move at threshold with identical velocities. The momenta of the baryons and the antiproton 
do not necessarily have to be equal because the effective masses of the nucleons and the 



24 



antiproton can differ from each other. In order to guarantee momentum conservation we 
use the following ansatz for the effective masses of the baryons and the antiprotons 

l''''l^ 3m'(|p,|)+.m; l'"+''^l- f'*" 

Since each baryon mass depends itself on the momentum of the corresponding baryon these 
equations are iterative equations for the absolute value of the momenta. Using these effective 
masses the threshold for antiproton production in the LRF reads 



yi^ = ^(3no + no + A0)2 - (sn,- + n^y. (57) 

p-production takes place if the invariant energy ^/s {s = {U^ + n^)^ — (11^ + n^)^) of the 
baryons in the initial channel lies above the threshold a/sq. According to definition (53) 



the elementary antiproton production cross section is a function of the invariant energy v s' 
available for the particles in the final state of the production reaction 



v^ = v/(n!] + no + ng + w^y - (p, + p, + p, + pp)^ (58) 

Using (|53| ) and (|54D this reads 



V^ = v/(n? + no - A0)2 - (p, + p^f. (59) 

This equation shows that in order to calculate the production probability one has to know 
the momenta of the particles in the final state of the elementary production reaction when 
employing momentum dependent potentials Us{x, p) and t/^(x, p) for the baryons. 

Since we are dealing with subthreshold particle production the main contribution to the 
production cross section will arise from events with collinear baryon momenta because this 
kinematical situation minimizes the kinetic energy of the baryons. Based on this argument 
we assume for the momenta of the baryons in the final state 

Pi = g(Pi+P2-Pp). (60) 
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We have applied the above mentioned formahsni to evaluate the antiproton cross section 
for the reactions p +^^ C and p+^^Cu at bombarding energies of 5, 4, and 3.5 GeV. The 



corresponding invariant cross sections in comparison with the data of ref. |Tl| are shown 
in fig. 18 as a function of the momentum of the emitted antiproton in the lab. frame at 
6 = 0°, assuming free antiprotons, i.e. g^, g^ = 0. The calculations slightly underestimate 
the experimental data, but already approximately reproduce the shape of the momentum- 
spectra as well as the dependence on bombarding energy and mass. The p reabsorption 
amounts to a factor of 12 in case of ^^C and to a factor of 19 for ^^Cm roughly in line with 
simple geometrical estimates. 

When adjusting the constant g^ such that the scalar potential (|T8|) becomes slightly 
attractive ( —50 to —100 MeV at po) the reproduction of the data improves at all ener- 
gies significantly which is exemplified for 4.0 GeV by the dashed line in fig. 18. In the 
above comparison we cannot distinguish between scalar and vector antiproton selfenergies 
because both yield similar results for the p spectrum if the same Schrodinger-equivalent 
optical potential is achieved. Furthermore, when using antiproton selfenergies in line with 
the relativistic mean- field theories [^, i.e. changing only the sign of the nucleon vector 
potential, we overestimate the p yield by more than an order of magnitude at all energies 
for both systems. 

We now turn to the nucleus-nucleus case. The calculated antiproton invariant differential 
cross section for the reaction "^^Si +^^ Si at 2.1 GeV/A and Ni -|- Ni at 1.85 GeV/u is 



shown in fig. 19 in comparison to the experimental data of ref. and ref. [0; the upper 



lines represent the results of the calculations for free antiprotons without including any 
reabsorption. When taking care of antiproton annihilation according to eq. ([A2|) the yields 



drop to the lower full lines. In the case of Ni + Ni the data are now underestimated sizeably. 
However, using attractive scalar (or vector) selfenergies at p = po of about —100 to —150 
MeV we reproduce the data for Ni + Ni, for Si + Si, however, we miss the data point at 
p = 1 GeV/c. 

We thus use a potential of —100 MeV at p = po to predict the differential p-excitation 
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function in Ni + Ni collisions from 1.4 to 2.5 GeV/u (fig. 20), a system that will be explored 



at GSI in the near future P2 |. 



The different value for the attractive antiproton field at p = po in p + A and A + A 
reactions is due to the fact that in p + A collisions the antiprotons move with momenta of 1 
- 2 GeV/c with respect to the nuclear medium, whereas in A + A collisions the antiprotons 
have smaller momenta in the nucleus-nucleus center-of-mass frame. In view of uncertainties 
of our present studies with respect to the elementary production cross sections close to the 
thresholds we provide areas for the antiproton Schrodinger equivalent potential at p = po 
in fig. 21, as extracted from the comparison with the experimental data for p + A and 
y4 + A reactions. These areas are far from the values expected by charge conjugation from 
the familiar a — u model [Q (dashed line) and thus exclude relativistic mean-field models 
with the same parameter-sets for nucleons and antinucleons. However, our extracted values 
are well in line with a Schrodinger-equivalent potential (solid line in fig. 21) as calculated 
from the dispersion relation (see equations (|1|) - (^). Crucial for this result is the correct 
momentum dependence for the p + A potential in the entrance channel. If we use the 
original Walecka model without the momentum dependent coupling strength we are forced 
to compensate for the strong repulsion with much deeper antiproton potentials. 

The primordial antiproton production rates and the reduction due to absorption obtained 
here agree well with those obtained by Li et al. p5|. These authors also get agreement 
with the data for Si + Si, however, by using the very deep p-potential obtained by charge 
conjugation from the Walecka model and the corresponding strong repulsion in the entrance 
channel; using the same model we also reproduce the data, as well as their calculations within 
a factor of two. However, the entrance channel potential is much too repulsive, as mentioned 
in section 3.1, so that an unphysically deep antiproton potential is needed to compensate for 
this repulsion. This fact shows that the conclusion of ref. [^ that the agreement of their 
calculation with the data supports the assumption of a very deep G-parity transformed 
nucleon potential for the antiprotons is not justified. 
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VII. SUMMARY 

In this work we have evaluated the differential cross section for p production for proton- 
nucleus and nucleus-nucleus reactions in the subtreshold regime by considering on-shell 
baryon-baryon production channels involving nucleons and A's with their in-medium quasi- 
particle properties and treating p propagation and annihilation nonperturbatively. The 
quasi-particle properties of the nucleons are fixed in our approach by the nuclear saturation 
properties, the proton-nucleus empirical potential as well as Dirac-Brueckner calculations 
at higher density. By varying especially the antiproton selfenergy we have shown that 
the differential p spectra are highly sensitive to the antiproton quasi-particle properties. 
Though we still have to rely on proper extrapolations to threshold of the elementary process 
p + p —>■ p + X, the latter sensitivity can be used to obtain approximate values for the p 
Schrodinger-equivalent potential from a systematic comparison to the available experimental 
data. 

In this respect we have performed a systematic study of p-nucleus and nucleus-nucleus 
collisions in a broad kinematical regime and compared our numerical results to the data 
from KEK [jll| and GSI [|12|. We find a consistent description of all data employing a 



rather weak attractive potential for the antiprotons which is well in line with a dispersive 
potential extracted from the dominant imaginary part of the antiproton selfenergy due to 
annihilation. Essential for this result is the use of the correct momentum dependence of 
the nucleon-nucleus potential in the entrance channel. Although the validity of the simple 
dispersive model used here may be questionable at higher densities and the extrapolation 
of the elementary p-production cross section down to threshold introduces an uncertainty, 
the consistency of the p-potentials obtained from a fit to the data with those calculated in 
the dispersive model indicates that the basic antiproton production mechanism in p + A and 
A + A reactions is understood. 

It has become clear that the baryonic production channels at subthreshold energies in- 
volve dominantly one or two nucleon resonances and that the production proceeds at the 
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highest densities that can be reached in a nucleus-nucleus reaction. Only at these high den- 
sities the relative population of nucleon resonances (~ 30 % as shown in |^) as well as the 
nucleon-resonance collision rate are high enough to allow for p production. Otherwise the 
resonances decay to a nucleon and a pion before colliding with another baryon such that the 
energy stored in the resonance gets lost for the production. 

We note in closing that the antiproton production studies at the AGS p^- ^] around 15 



GeV/u, although far above the free production threshold, might yield further information 
on the dynamics and selfenergies of antiprotons at even higher densities. 

The authors acknowlegde valuable discussions with A. Gillitzer, P. Kienle, W. Konig and 
A. Schroter as well as their information on experimental results prior to publication. They 
are also very grateful to C. M. Ko for detailed discussions both on the physics and on the 
computational aspects of these calculations. 
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APPENDIX A: THE WEIGHTED TESTPARTICLE METHOD 

In the following we present the weighted testparticle method used to implement the 
antiproton production, absorption and propagation in the RBUU model. It is especially 
suited for the simulation of processes with low production probabilities and/or high in- 
medium absorption cross sections for which conventional testparticle methods fail because 
of low statistics. 

Assuming that during a calculation N elementary baryon-baryon collisions with sufficient 
energy for antiproton production occur, we define for each coUison i (ie{l, ■■■, ,N}) a 
three-dimensional momentum grid in the reference frame (CMS of the heavy-ion collision 
or lab-system). In order to calculate the production probability for each collision i we 
transform this grid into the system that fuUfiUs eq. (|4^) or into the LRF (see eqn. (|49|) to 
(|60|)). After calculating the invariant differential p-production probability wj!/^ for each grid 
point p= {kApx, lApy, mApz) the momentum grid is transformed back into the reference 
frame. To obtain the total invariant differential production probability without absorption 
for a heavy-ion collision or a proton-nucleus reaction one has to sum up the contributions 
from all the N elementary production reactions 

N 



a 

i=l 



WMm = Y.<L (Ai) 



In order to calculate the antiproton probability including propagation through the dense 
medium and absoption we start out with the momentum grid described above. For all N 
baryon-baryon collisions contributing to the antiproton production we place a testparticle 
representing an antiproton on every point of the momentum grid. This yields additional j 
{j = k ■ I ■ m) testparticles for each of the N baryon-baryon collison events. The effective 
masses and momenta of these testparticles are determined according to both the location 
on the momentum grid and the spatial coordinate where the baryon-baryon collision takes 
place. In addition we attach variable weights to the testparticles. The weight W;!' of the 
h-th testparticles 'created' in the i-th collision denotes the invariant differential production 
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probability tf^/m ^^ ^^^ corresponding grid point. The result of this procedure is an ensemble 
of testparticles representing the produced antiprotons that contain information about the 
differential p-distribution for each relevant baryon-baryon collision. The entire testparticle 
distribution together with the corresponding weights represents the phase-space distribution 
of the antiprotons produced in a heavy-ion collision or a proton-nucleus reaction. 

We now describe how we treat this ensemble of testparticles in order to calculate the 
antiproton propagation and absorption: Each testparticle is propagated by means of the 
equation of motion for the antiparticle phase-space distribution function (^). At the end 
of each timestep the absorption rate for the testparticles is calculated according to their 
weights. In order to perform this task we convert the collision integral (|33|) into a differential 
equation for the testparticle weights 



dt ' ' " (2 
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5\W^ + n^ - n:^)/(a:^ n^)«^f U,n„ (A2) 

where x^ is the space-time coordinate where the annihilation takes place. lip, II^, nip and 
m*B denote the effective momenta and masses of the antiprotons and the baryons, respec- 
tively. This differential equation leads to the following changes in the testparticle weights 
in timesteps At, 

= w'{t,)fexp{- f'"' j ^m^^W{p + B ^ X)f{x^,Ii%)dt). (A3) 

Solving the integral in the exponent of (^5) by means of the local density approximation 
(IDA) leads to 

nj\t,^i)f = ^'(^.)?e^P(-^ E I § - ^ I ^iP + Bj - X)At), (A4) 

with Nt being the number of testparticles per nucleon used to model the baryon phase-space 
distribution and AV being the normalization volume used in the RBUU model to evaluate 
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densities and currents. Eq. (^J) allows for a nonperturbative evaluation of high absorption 
rates without a reduction of the number of testparticles. 

At the end of the simulation we project the weight of each antiproton testparticle obtained 
in the final timestep tf of our calculation on the three-dimensional momentum grid (see 
above) 

w'h {tf, p^, py, p,) -^ w'li^ (A5) 

and sum up all contributions in order to obtain the invariant differential antiproton proba- 
bility for the corresponding heavy-ion collision or proton-nucleus reaction . 
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APPENDIX B: FIGURE CAPTIONS 



fig.l: Free annihilation cross section ioi p + p ^ X. Dots: experimental data from 



2S]; solid line: parametrization 



fig. 2: Real part of the antiproton selfenergy as a function of its kinetic energy using 
different parameters B in (^) for the annihilation cross section. 

fig.3: Inclusive antiproton cross section in a pp collision. The solid line indicates the 
parametrization (^7\). The dashed and dotted lines show 'extreme' alternative 



extrapolations to the threshold. The experimental data are taken from [T^ 



fig.4: Invariant differential production probability for Si + Si at 2.1 GeV/u in the CMS 
(e = 0°; b = 0; EOS: K = 200 MeV, m*/m = 0.7 at p = po) using different 
parametrizations for the elementary production cross section according to fig. 3. 

fig.5 : Effects of Pauli-blocking for Au + Au at 2.1 GeV/u in the CMS (6 = 0°, 6 = 0); 
Solid line: Pauli-blocking included; dashed line: no Pauli-blocking. 

fig.6: Differential p-production probability for Si + Si (top) andA-u + Au (bottom) at 
2.1 GeV/u in the CMS (O = 0", 6 = 0) using different parametrizations for the 
EOS; NLl: K = 400 MeV, m* jm = 0.83; NL2: K = 200 MeV, m*/m = 0.83; 
NL3: K = 200 MeV, m*/m = 0.83 at p = po- 

fig. 7: Production probability as a function of p/po for Au + Au (solid line) and Si + Si 
(dashed line; enhanced by a factor of 50) for 2.1 GeV/u; Qcm = 0°, b = 0. 

fig.8: p-multiplicity for Si + Si at 2.1 GeV/u as a function of the impact parameter 
b weighted with 27r6; calculation: solid line; analytic formula: dashed line (cf. 

text). 

fig.9: Contributions to the production from different channels: NN (filed histograms), 
A^A (open histograms) and A A (hatched histograms) for Si + Si, Ca + Ca and 
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Au + Au at 2.0 GeV/u and 2.5 GeV/u. 

fig. 10: Relative contributions to the p-production from different cliannels as a function 
of tlie beam energy for Ca + Ca; NN (dotted line), A^A (solid line) and A A 
(dashed line). 

fig. 11: Differential p-production probability for a central Ni + Ni collision at 1.85 
GeV/u {QcMS = 0°) as a function of the antiproton momentum. 
Top: fixed U^ = -98 MeV at p = po and variing UP at p = po: -159 MeV 
(solid), -100 MeV (dashed), -40 MeV (dotted) and (dashed-dotted). 
Bottom: fixed U^ = -159 MeV at p = po and variing U^ at p = po: -98 MeV 
(solid), -72 MeV (dashed), (dotted) and +72 MeV (dashed-dotted). 

fig. 12: Effects of the in-medium propagation on the p-sprectrum oi Si + Si at 2 A 
GeV/u under 0° in the CMS. 
Top: UP = and U^ = -72 MeV (dashed), (solid) and +72 MeV (dotted) at 

P = Po- 

Bottom: f/| = and t/f = -100 MeV (dashed), (sohd) and +100 MeV 

(dotted) at p = po. 

fig. 13: Effects of the in-medium propagation on the p-spectrum for the reaction p + Cm 
at 4.0 GeV, b= 0, Qiab = 0° for different optical potentials: Usep = (dotted), 
Usep = -90 MeV (solid) and Usep = -560 MeV (dashed). 

fig. 14: Number of antiproton -baryon collisions for a central Si + Si reaction at 2.1 
GeV/u as a function of the corresponding elementary p-annihilation cross section. 

fig. 15: Antiproton spectra for the reaction Si + Si at 2.1 GeV/u (b = 0, Qcms = 0°) 
using different elementary annihilation cross sections: upper line: no absorp- 
tion; dotted line: aabs = 60 mb; dashed-dotted line: aabs = 70 mb; solid line: 
parametrization (^; dashed line: aabs = 80 mb. 
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fig. 16: Quotient R as a function of the impact parameter scaled with bmax for the 
systems Ni + Ni, Ne + Ni and Si + Si at 1.85 GeV/u, 2.0 GeV/u and 2.1 
GeV/u. Dots: RBUU- model; solid lines: fits according to eq. (^). 

fig. 17: Scalar and vector selfenergies Us{p) and Uq{p) for nucleons at different densities 
in units of po ~ 0.17 fm~^. 

fig. 18: Invariant cross section for antiproton production in the reactions p+^^C and 
P_I_63q^ at O = 0° as a function of the antiproton momentum p in the lab. 
system. The experimental data are taken from ref. [ill] and correspond to bom- 



barding energies of 5.0 GeV , 4.0 GeV and 3.5 GeV. The full lines represent 
calculations for free antiprotons. The dashed lines indicate the result for an 
antiproton selfenergy of - 100 MeV at 4.0 GeV. 

fig. 19: Invariant cross section for antiproton production in the reaction '^^Si +^^ Si at 
2.1 GeV/u and Ni + Ni at 1.85 GeV/u for = 0° as a function of the momentum 
of the emitted antiproton in the lab-system. The experimental data have been 
taken from refs. [Q and |[T^ , respectively. The upper lines indicate the calculated 
cross section for free antiprotons without reabsorption whereas the lower solid 
line is obtained when including p annihilation. The dashed line represents the 
cross section adopting an attractive potential of the antiproton of - 150 MeV. 

fig. 20: Invariant cross section for antiproton production in the reaction Ni + Ni at 2.5 
GeV/u (sohd), 2.0 GeV/u (dashed), 1.85 GeV/u (dotted), 1.6 GeV/u (dotted- 
dashed) and 1.4 GeV/u (dotted-dotted-dashed). All calculations were done using 
an attractive p-potential of -100 MeV at p = po. The experimental data are taken 
from iT^ . 



fig. 21: Comparison of our extracted values for the Schrodinger equivalent antiproton 
potential from p + A and A -|- A reactions with the prediction from the a — u 
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model (dashed line) and the dispersive potential according to eqns. (||) - (|) 
(solid line). 
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TABLES 



2.0 GeV/u 



2.5 GeV/u 



NN 



NA 



AA 



NN 



NA 



AA 



Si+Si 


5 % 


45 % 


50 % 


Ca+Ca 


6 % 


46 % 


48 % 


Au+Au 


4 % 


42 % 


54 % 



10 % 



6 % 



52 % 


38% 


50 % 


42 % 


44 % 


50 % 



TABLE L Relative contributions of the different reaction channels integrated over Ni + ^"2- 



Ni + Ni 



Ne + Ni 



Si + Si 



A 
B 
C 



1.0022 

0.9 

0.05 



1.0025 

0.9 

0.04 



1.0122 

0.9 

0.06 



TABLE IL Parameters A, B and C for the fit m) of R. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9401008vl 



This figure "fig2-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9401008vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9401008vl 



This figure "fig2-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9401008vl 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9401008vl 



